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REVIEW  OVER  COMPLETE  GRANT  PERIOD  1972-1976 


1.1  Major  Developments 

Most  of  the  research  has  been  associated  with  the  further 
development  and  application  of  the  direct  simulation  Monte  Carlo  method 
for  the  modeling  of  gas  flows  at  the  molecular  level.  A monograph1  which 
includes  a full  description  of  the  method,  together  with  a number  of  new 
applications,  has  been  written  and  published  during  this  period.  Acknow- 
ledgment was  made  generally  of  the  support  provided  by  the  AFOSR  during  the 
whole  period  over  which  the  method  has  been  developed  and  specifically  to 
the  current  grant  with  respect  to  some  new  applications  that  appear  in 
the  book. 

A significant  development  has  been  the  demonstration  that 
the  computations  are  best  carried  out  on  a dedicated  minicomputer  rather 
than  on  a large  machine  in  a computing  center.  The  probability  of  this 
was  proposed  for  preliminary  investigation  in  1973  and  this  led  to  a 
proposal  to  the  Australian  Research  Grants  Comnittee  for  a PDP  11/40 
system.  This  submission  was  successful  and  the  computer  was  installed 
early  in  1975.  Since  this  time,  it  has  run  twenty  four  hours  a day  and 
seven  days  a week  without  an  operator  and  the  improvement  in  cost-effect- 
iveness has  been  found  to  be  near  two  orders  of  magnitude.  The  reduction 
in  computing  costs  has  meant  that  small  disturbance  flows,  which  had 
previously  been  thought  to  be  inaccessible  to  the  method,  can  now  be 
studied.  The  first  demonstration  of  this  capability  was  for  the  classical 
surface  slip  or  Kramers  problem*.  The  flow  profile  was  studied  at  Mach 
numbers  as  low  as  0.05,  with  the  computer  runs  being  sufficiently  long 
to  reduce  the  statistical  scatter  to  one  part  in  a thousand. 

An  additional  consequence  of  the  practicability  of  very  long 
runs  is  that  molecular  simulation  can  now  be  considered  is  an  alternative 
to  finite  difference  methods  for  the  computation  of  some  continuum  flows. 

This  point  was  first  demonstrated4  by  lengthy  runs  on  a large  computer  and 
a more  recent  study*  has  further  demonstrated  the  advantages  of  the  dedicated 
minicomputer.  The  advantages  of  the  molecular  simulation  over  the  finite 
difference  continuum  representation  relate  to  the  absence  of  stability 
problems  and  the  ease  of  incorporation  of  conplex  flow  boundaries.  The 
molecular  approach  appears  to  be  most  advantageous  for  problems  involving 
complex  shock  wave  interactions,  but  having  inviscid  boundaries.  On  the 
other  hand,  the  method  has  been  used5  to  study  a viscous  problem  at 
sufficiently  low  Knudsen  number  that  the  conclusions  should -be  valid  in 
continuum  flow. 

The  most  recent  extension  of  the  method  has  been  to  reacting 
gas  flows.  Particular  attention  has  been  paid  to  the  dissociation- 
recombination  reaction  in  a diatomic  gas.  The  model  permits  the  analytical 
determination  of  the  equilibrium  state1  and  it  has  now  been  applied  to  the 
reaction  zone  behind  a strong  shock  wave.  Specific  cases  have  been  computed 


for  nitrogen  and  the  results  have  been  compared  with  existing  experimental 
data  for  nitrogen10. 


One  topic,  the  study  of  the  forces  on  very  small  suspended 
particles  in  the  atmosphere  started*  as  a simulation  project,  but  led  to 
a significant  experimental  study6. 
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2,  REPORT  ON  PERIOD  1975-1976 


The  final  extension  of  the  Grant  was  from  April  15,  1975 
to  April  14,  1976  but,  due  to  a period  of  absence  of  the  Principal 
Investigator,  this  was  extended  to  October  14,  1976.  Four  topics  were 
included  in  the  proposal  and  the  work  done  on  each  of  these  is  described 
below. 


2.1  Collision  Rates  and  Colllsional  Energy  Distribution  within 
Shock  Waves. 


This  work  lad  to  publication  ? in  the  list  of  Section  1.2 
The  direct  simulation  Monte  Carlo  method  was  used  to  determine  the 
collision  rate  and  the  distribution  of  relative  translational  energy 
between  colliding  molecules  within  strong  shock  waves.  The  results 
were  compared  with  the  corresponding  theoretical  quantities  in  an 
equilibrium  gas.  The  objective  was  to  determine  whether  the  region  of 
marked  translational  nonequilibrium  within  the  wave  needs  to  be  taken 
into  account  in  the  interpretation  of  reaction  rate  experiments  in 
shock  tubes.  It  was  found  that  a correction  should  be  applied  to 
measured  relaxation  distances  behind  shock  waves.  The  magnitude  of 
this  correction  is  of  the  order  of  ten  times  the  translational  shock 
thickness. 


2.2  Calculations  with  Quantum  Mechanics  Cross-sections 


Dr.  A.U.  Chatwani  of  the  Max-Planck  Institut  fur 
Stromungsforschung  spent  six  weeks  in  the  Aeronautical  Engineering 
Department  and  incorporated  quantum  theory  collision  cross-sections 
into  the  existing  simulation  program  for  a steady  spherical  expansion. 
This  program  had  been  used  to  study  the  onset  of  translational  non- 
equilibrium, and  this  generally  occurs  at  very  low  temperatures  where 
quantum  effects  are  significant.  Dr.  Chatwani  continued  the  work  after 
his  return  to  Germany  and  carried  out  a thorough  study  of  terminal  Mach 
numbers,  both  for  real  gases  and  the  idealised  models.  This  work  was 
presented  at  the  Tenth  International  Symposium  on  Rarefied  Gas  Dynamics 
in  a paper  entitled  "Monte  Carlo  Simulation  of  Nozzle  Beam  Expansions" 
and  is  to  be  published  in  the  proceedings. 


2.3  The  Onset  of  Mach  Reflection 


The  direct  simulation  Monte  Carlo  method  was  used  to  study 
the  transition  from  regular  to  Mach  reflection  of  oblique  two-dimensional 
steady  shock  waves.  Particular  attention  was  paid  to  the  region  of  Mach 
number  and  deflection  angle  in  which  both  types  of  reflection  are 
theoretically  possible.  The  results  indicated  that  regular  reflection 
can  occur  in  the  region  where  both  types  of  reflection  are  possible.  The 
study  confirmed  that  the  molecular  approach  can  provide  a practical 
alternative  to  finite  difference  and  finite  element  methods  for  the 
computation  of  some  continuum  flows,  especially  when  a dedicated  mini- 
conputer  is  used  instead  of  a large  machine  in  a general  coaputing  center. 
A paper  was  presented  at  the  AIAA  9th  Fluid  and  Plasma  Dynamics  Conference 


2.4  The  Dissociation-Recombination  Reaction  in  a Shock 


The  direct  simulation  Monte  Carlo  method  has  been  extended 
to  cover  molecular  dissociation  and  recombination  reactions.  The  model 
also  includes  internal  degrees  of  freedom  with  separate  relaxation  times 
for  rotation  and  vibration.  The  model  is  based  on  a classical  collision 
theory  analysis  that  provides  an  expression  for  the  equilibrium  state  and 
allows  the  appropriate  simulation  parameters  to  be  inferred  from  thermo- 
dynamic data  for  real  gases.  The  model  was  first  tested  by  the  computation 
of  the  constant  volume  relaxation  of  a homogeneous  gas  sample.  A finite 
difference  procedure  was  then  used  to  extend  this  program  to  cover  the 
dissociation  relaxation  zone  behind  strong  normal  shock  waves.  A number 
of  nitrogen  cases  were  computed  and  the  results  were  in  agreement  with 
existing  experimental  data.  A paper  has  been  submitted  to  the  Journal  of 
Computational  Physics. 


In  addition  to  the  work  on  the  above  topics,  the  work  on 
the  incompressible  Kramers  problem  (reference  s in  the  list  of  Section 
1.2)  was  carried  out  during  this  period. 
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